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ABSTRACT 

We report the identification of the M9 dwarf SDSS J000649.16-085246.3 as a spectral binary and 
radial velocity variable with components straddling the hydrogen burning mass limit. Low-resolution 
near-infrared spectroscopy reveals spectral features indicative of a T dwarf companion, and spectral 
template fitting yields component types of M8.5±0.5 and T5±l. High-resolution near-infrared spec- 
troscopy with Keck/NIRSPEC reveals pronounced radial velocity variations with a semi-amplitude 
of 8.2±0.4 km s _1 . From these we determine an orbital period of 147.6±1.5 days and eccentricity 
of 0.10±0.07, making SDSS J0006-0852AB the third tightest very low mass binary known. This 
system is also found to have a common proper motion companion, the inactive M7 dwarf LP 704-48, 
at a projected separation of 820±120 AU. The lack of Ha emission in both M dwarf components 
indicates that this system is relatively old, as confirmed by evolutionary model analysis of the tight 
binary. LP 704-48/SDSS J0006-0852AB is the lowest-mass confirmed triple identified to date, and 
one of only seven candidate and confirmed triples with total masses below 0.3 M Q currently known. 
We show that current star and brown dwarf formation models cannot produce triple systems like 
LP 704-48/SDSS J0006— 0852AB, and we rule out Kozai-Lidov perturbations and tidal circulariza- 
tion as a viable mechanism to shrink the inner orbit. The similarities between this system and the 
recently uncovered low-mass eclipsing triples NLTT 41135AB/41136 and LHS 6343ABC suggest that 
substcllar tertiaries may be common in wide VI dwarf pairs. 

Subject headings: binaries: spectroscopic — binaries: visual — stars: fundamental parameters — 
stars: individual (SDSS J000649.16-085246.3, LP 704-48) — stars: low mass, 
brown dwarfs 



1. INTRODUCTION 

Our current theoretical and observational understand- 
ing of star formation holds that most stars form in mul- 
tiple systems, with the frequency, mass ratio distribu- 
tion and period distribution of multiples varying as a 
functi on of mass and possibly formation environment 
(e.g.. IDelgado-Donate et all 12004 lAhmic et all 120071: 
Burgasser et all l2007bl: iBatd 120091 120121 iKraus et al.l 



20111 ). Dynamics play an important role in multi- 



ple formation, as both the fragmentation of collaps- 
ing cl oud cores (IBossI 120011) and massive circumstellar 
disks (jStamatellos fc Whitworthll2009t IStamatellos et al.l 
1201 If ) give rise to small N groups t hat either dissolve 
or sta bilize into hierarchical systems ()Sterzik fc Durisenl 
2003). Numerical calculations show that both fragmen- 
tation and dynamical evolution of stellar multiples are 
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sensitive to the initial cloud properties — geometry, den- 
sity, turbulence spectrum, and radiative and external 
feedback — so multiplicity statistics provide important 
tests for star formation theory, on par with the initial 
mass function, stellar mass segregation and velocity dis- 
tributions. 

The mass dependence of multiplicity among stars 
is well established, with observ ed frequencies rang- 
ing f ro m >80% for OBA dwarfs (IS hat sky fe Tokovinin 
20021: iKouwenhoven et all 120051: iKobulnickv fc Fryer 
30% for sy s tems with M dwa r f pri- 
Marcvl [19921 : IReid fc Gizisl fl997t 
The frequency of multiples ap- 
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maries dFischer _. 
IDelfosse et al.l 120041) 
pears to drop even further in the very low mass (VLM; 
M < 0.1 M©) stellar and substellar regimes. Re- 
solved imaging studies of field and cluster VLM dwarfs, 
based on ground-based adaptive optics (AO) and Hub- 
ble Space Telescope (HST) observations, find relatively 
consistent (volume-limited) binary frequencies of 10- 
20% Ce.g. lBouv et al.ll20MlBurgasser et al.ll200l l2006d : 
iClose et al.ll2003l) . However, estimates for the frequency 
of tightly-bound spectroscopic VLM binaries (roughly 
10% of currently known systems ) span a broad range of 
1-25% (iMaxted fc Jeffries] 120051: iBasri fc Reinersl |200 
Kurosawa et al.ll2006t I Joergensl I2OO8I: iBlake et aLll201 
Clark et al.l 120121 ) . As such, there remains significant 



uncertainty in the overall VLM binary frequency, which 
has a direct impact on our under standing of how brown 
dwar fs form in the first place (e.g.. IDelgado-Donate et all 
2003). There is also evidence that the orbital characteris- 
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tics of multiples are mass-dependent, with VLM binaries 
being on average tighter ((a) w 7 AU versus w 30 AU) 
and more frequently composed of equal-mass components 
((q) = M2 /M i » 1 versus ~ 0.3) than their solar-mass 
counterparts (Allen] I2007D . Again, these trends may 
be skewed by biases inherent to the known sample of 
VLM binaries, identified largely as imaged pairs. A com- 
plete understanding of VLM multiplicity statistics there- 
fore requires more robust constraints on the frequency 
and characteristics of s hort-period and \ow-q multiples 
(jBurgasser et al.ll2007bh . 

An alternative method for identifying VLM multiples 
that avoids separation biases (both intrinsic and pro- 
jected) is through spectral binaries, systems with com- 
ponents of different spectral types whose combined-light 
spectra exhibit distinct peculiarities. While this tech- 
nique is more comm only associated with white dwarf- 
M dwarf pairs (e.g., iSilvestri et al.l 120061 ) . spectral bi- 
naries containing late M or L dwarf primaries and T 
dwarf secondaries have also been recognized due to 
thei r unique and highly structured near-infrared spec- 
tra dCruz et all 120041: iBurgasserl 1200711 iBurgasser et all 
2008al iStumpf et al l l2008t JBurgasser et al.ll2010L I2011ri 
Gelino fc BurgasserlboioT iGeiffler et , all 1201 lft . Because 
blended light systems can be identified at any separa- 
tion, M/L+T spectral binaries can probe the very clos- 
est separations that are inaccessible to direct imaging 
studies (a < 50 mas). An illustrative case is the M8.5 
+ T5 binary 2MASS J03202839-0446358 (hereafter 
2MASS J0320-0446), a sys tem independently iden tified 
as both a spectral binary (Burgass er et all [2008aD and 
a radial velocity (RV) variable with a period of 0.68 yr 
and primary semi- major axis a\ sini = 0.157±0.003 AU 
([Blake et all 120081 12010ft . With a projected separation 
of «17 mas, this system is unresolvable with current 
imaging and interferometric technology. Moreover, the 
significant difference in component magnitudes (AK = 
4.3±0.6) makes this system a potential low q pair. The 
independent constraints provided by the component clas- 
sifications and RV orbit yield robust limits on the age 
(>2 Gyr) and orbit al inclination (> 53°) of this system 
([Burgasser fc Blakd [2009ft . 

In an effort to confirm and characterize spec- 
tral binaries identified in low-resolution near- 
infrared spectroscopy, we have identified a new 
VLM system exhibiting significant RV variability. 
The source, SPSS J000649. 16-08 5246.3 (hereafter 
SDSS J0006-0852; IWest et al.l [2008ft is an optically- 
classified M9 which appears to be both a short-period 
(0.4 yr) binary with M8.5 and T5±l components, 
and a co-moving wide companion to the inactive M7 
dwarf LP 704-48. Together, this system comprises the 
lowest-mass triple confirmed to date. In Section 2 we 
describe our optical and near-infrared imaging and 
spectroscopic observations of this system using the lm 
Lick Observatory Nickel Direct Imaging Camera, the 
3m NASA In frared Telescope F acility (IRTF) SpeX 
spectrograph (Rayncr et all 12003ft . the 4m Kitt Peak 
National Observatory (KPNO) RC spec trograph, and 
the 1 0m Keck II NIRSPEC spectrograph ([McLean et all 
1998). In Section 3 we examine the properties of the 
SDSS J0006-0852 and LP 704-48 pair, assessing their 
common proper motion, distance and probability of 
chance alignment; as well as the age and metallicity 
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Fig. 1 — Lick/Nickel image of the LP 704-48/SDSS J0006-0852 
field showing the location of the two sources (arrows and labels) and 
field sources used for calibrating the astrometric reference frame 
(circles). The field shown is approximately 5'x5' and oriented 
North up and East to the left. Streaking along the bottom of the 
image is due to reflected lunar light. 

of both components as derived from spectroscopic and 
kinematic indicators. In Section 4 we analyze the low- 
resolution near-infrared spectrum of SDSS J0006— 0852 
to infer the presence of, and characterize, its T dwarf 
companion. In Section 5 we analyze our RV mea- 
surements which allow us to extract both the orbital 
properties (including constraints on the orbital inclina- 
tion) and verify the relatively old age for this system. In 
Section 6 we discuss LP 704-48/SDSS J0006-0852AB 
in the context of other hierarchical low-mass triples, 
and examine whether current star formation theories or 
three-body interactions can create such systems. Our 
results are summarized in Section 7. 

2. OBSERVATIONS 
2.1. Optical Imaging 

The SDSS J0006-0852 field was observed with the 
Nickel lm CCD camera on 2011 October 28 (UT) in 
clear conditions with l'/5 seeing. The chip was binned 
2x2 for a pixel scale of 0'.'37 pixel -1 . Six dithered ex- 
posures of 150 s each were obtained using the /-band 
filter over an airmass range of 1.45-1.49. Dome flat and 
dark exposures were obtained at the end of the night for 
pixel calibration. Data were reduced using custom IDlQ 
scripts that removed the bias voltage, median-combined 
and normalized the flat field frames; divided the science 
data by the flat; masked bad pixels; generated a median 
sky frame that was subtracted from the science data; 
and registered and median-stacked the science frames to 
a single final image, shown in Figure [TJ 

Pixel positions for sources in the final stacked frame 
were measured using a cent roiding algor ithm routine at- 
tached to the atv program (Barthl [2001ft . Using astrom- 
etry of 13 field stars drawn from the Two Micron All 
Sky Survey (2MASS; iSkrutskie et~aTI 12006ft . we deter- 
mined a linear astrometric transformation frame with 
residuals of 0'.'3 in Right Ascension and 0'.'2 and decli- 
nation; i.e., ^20% of the seeing radius. This transforma- 
tion was used to compute coordinates for LP 704-48 and 
SDSS J0006-0852, listed in Table LU Figure M displays 
the 2MASS to Nickel offsets for all of the sources in the 

7 Interactive Data Language. 
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Fig. 2. — Offsets in Right Ascension and declination between 
2MASS and Nickel imaging epochs for field stars (plus symbols) 
and LP 704-48/SDSS J0006-0852 (diamonds). The error bars on 
the last two data points reflect the scatter in offsets among the 
field stars. 



TABLE 1 

Astrometry for SDSS J0006-0852 and LP 704-48 



Source 


Epoch 


a 


8 


Uncertainties 




CUT) 


(J2000) 


(J2000) 


(") 




LP 704-48 


POSS I R 


1954.669 


00 06 47.66 


-08 52 21.17 


0.1, 0.1" 


2MASS 


1998.767 


00 06 47.46 


-08 52 35.00 


0.11, 0.09 


ESO R 


1999.604 


00 06 47.45 


-08 52 35.31 


0.1, 0.1 a 


UKST I N 


2000.723 


00 06 47.45 


-08 52 35.66 


0.1, 0.1 a 


SDSS 


2000.740 


00 06 47.45 


-08 52 35.60 


0.06, 0.04 


Nickel 


2011.820 


00 06 47.39 


-08 52 39.16 


0.3, 0.2 






SD9S J0006- 


-0852 




2MASS 


1998.767 


00 06 49.16 


-08 52 45.70 


0.11, 0.09 


ESO R 


1999.604 


00 06 49.15 


-08 52 45.53 


0.1, 0.1 a 


UKST I N 


2000.723 


00 06 49.14 


-08 52 46.42 


0.1, 0.1 a 


SDSS 


2000.740 


00 06 49.16 


-08 52 46.30 


0.07, 0.06 


Nickel 


2011.820 


00 06 49.09 


-08 52 49.96 


0.3, 0.2 


References 


- POSS 


1 R, ESO R 


and UKST I N 


astrometry are 



from the SuperCosmos Sky Survey IHamblv e~a l. 2001c b a); 2MASS 
astr ometry are from the All-Sky Data Release Point Source Cata- 
log (Skrutskic ct al. 2006); SDSS astrometry are from Data Release 7 
(Abazaiian ct al. 2009). 

a Estimated. 

field. The large and similar proper motions for LP 704- 
48 and SDSS J0006-0852 are evident, distinct from the 
field by > 20a in declination. 

2.2. Low Resolution Optical Spectroscopy 

The optical spectrum of SDSS J 0006-0852 from the 
Sloan Digital Sky Survey (SDSS; lYork et all l2000h is 
shown in Figure El compared to an M9 SPSS spec- 
tral template froni LBochariski "etaLl (I2007bl) . The sim- 
ilarity _of_fchese spectra confirms the M9 classification 
from IWest et all ([2008D . although we note a slightly 
flatter 7500 A VO band and weaker Na I lines at 
8200 A for SDSS J0006-0852. Neither Li I absorp- 
tion (EW < 0.3 A) nor Ha emission (/ flo < 4 x 
10~ 18 erg s _1 cm -2 ) are detected in these data. Up- 
per limits on log 10 Lna/Lboi were calculated by finding 
the Ha equivalent width corres ponding to the lcr flux 
uncertainty, using th e x factor (jWalkowicz et al.l I2004t 
IWest fe Hawlevll2008D to convert from EW to L H a/L bo i 
and to propagate uncertainties. We find \og 1Q Lna/Lboi 
< -5.7 for SDSS J0006-0852. 

The nearby proper motion star LP 704-48 (aka 




Fig. 3. — Red optical spectra of LP 704-48 (top) and 
SDSS J0006-0852 (bottom, fr om SDSS) compared to M 7 and M9 
SDSS spectral templates from iBochanski ct al. (2007b). All spec- 
tra shown are smoothed to a common resolution of A/ A A = 1200. 
Data for the two dwarfs are scaled to their apparent SDSS i magni- 
tudes, while the templates are scaled to overlap in the 7400-7500 A 
(LP 704-48) and 8550-8600 A (SDSS J0006-0852) regions. Pri- 
mary atomic and molecular absorption features are labeled. The 
inset boxes shows close-ups of the 6520—6730 A region, showing the 
absence of both Ha emission and Li I absorption in both sources. 

2MAS S J00064746- 0852350;lLuyt^[198g) was classified 
M6 by I Cruz fc Reidl (|2002[ ) based on low-resolution opti- 
cal spectroscopy. We obtained a new optical spectrum of 
LP 704-48 on 2011 Nov 28 (UT) using the RC spectro- 
graph on the KPNO 4m telescope. Conditions were good, 
with seeing ~1" and excellent transparency. We used a 
l'/5 slit aligned to the parallactic angle, a 600 lines mm -1 
grating blazed at 7500 A (BL 420), and an order-blocking 
filter (OG 530) to obtain 5800-9100 A spectra at a reso- 
lution of ^2000 and dispersion of 1.74 A pixel -1 . A sin- 
gle 1200 s exposure was obtained at an airmass of 1.32. 
We observed the sdO flux standard HZ 4 for flux cali- 
bration, and obtained a set of quartz flats, bias frames, 
and arclamp exposures for pixel response and wavelength 
calibration. Data were reduced using standard routines 
in IRArfl 

The reduced spectrum for LP 704-48 is shown in 
Figure [3 compared to an M7 SDSS template from 
Bochanski et al.l (|2007bT) . Both our data and those of 
Cruz fc Reid ( 200J), most closely match this compari- 
son source, with the exception of a downturn in flux be- 
yond 8100 A caused by poor flux calibration in the red 
and strong fringing. An M7 classification is also gener- 
ally consistent with spectral index-based classifications 
us ing the VOl VP 2 Ti07 and Color-M indices defined 
in lLepine et ail Tabled}. Like SDSS J0006-0852, 

LP 704-48 exhibits no evidence of Li I absorption (EW 
< 0.13) or Ha emission (fna < 3 x 10~ 17 erg s -1 cm~ 2 ; 
login L-fjn/Lhni < — 6.2) in either our spectrum or that of 
I Cruz fc Reidl (|2002l ). As discussed in Section 3, the lack 
of nonthermal emission is consistent with a relatively old 
age for both sources. 

8 Image Reduction and Analysis Facility (Todv 1986). 



4 



TABLE 2 

Optical Spectral Indices for LP 704-48 and SDSS J0006-0852 



Index LP 704-48 LP 704-48 a SDSS J0006-0852 Ref 



Ti05 


0.208±0.005 




0.213±0.010 




0.33±0.03 




1 


CaH2 


0.256±0.004 




0.255±0.007 




0.35±0.03 




1 


CaH3 


0.558±0.007 




0.569±0.013 




0.61±0.03 




1 


VOl 


0.836±0.004 


(M6.5) 


0.826±0.008 


(M7) 


0.768±0.014 


(M9) 


2.3 


V02 


0.582±0.003 


(M6.5) 


0.610±0.006 


(M6) 


0.340±0.007 


(M9) 


3 


Ti07 


0.707±0.006 


(M6) 


0.723±0.009 


(M6) 


0.504±0.014 


(M8) 


3 


Color-M 


4.85±0.03 


(M7) 


5.23±0.06 


(M7) 


9.41±0.15 


(M9.5) 


3 


c 


1.002±0.010 




1.004±0.014 




0.96±0.06 




1 



References. — (ll jReid et al.l 11993) : (2) IHawlev et al.1 (200% ; (3) ILepine et al.l (120031) : 
(4) ILepine et a).l (1200711 . 

a Data from lCruz fc Reidl (120021) . 

Classifications based on the relations quantified in lLepine et al.l d2003l ). 



2.3. Low Resolution Near-Infrared Spectroscopy 

Low resolution near- infrared spectra of 
SDSS J0006-0852 were obtained on 2009 Novem- 
ber 4 (UT) with IRTF/SpeX in conditions of light cirrus 
and 0'.'6 seeing. We used the prism-dispersed mode 
with the 0'.'5 slit to obtain continuous, low-resolution 
(A/AA fa 120) spectra covering the 0.7-2.4 /mi band. 
Six exposures of 120 s each were acquired in an ABBA 
dither pattern along the slit which was aligned with 
the parallactic angle. The AO V star HD 1154 (V = 
8.85) was observed immediately afterward at a similar 
airmass for flux calibration and telluric absorption cor- 
rection, and NeAr arc lamp and flat field exposures were 
obtained along with the standard. Data were reduced 
with the SpeXtool pack age, version 3.4 (dishin g et al.l 
12004 iVacca et ail I2003T) . using standard settings. A 
de tailed desc r iption of our reduction procedures is given 
in iBurgasserl (|2007bh . 

Figure H] displays the reduced spectrum of 
SDSS J0006— 0852 compared to equi valent data for 
2MAS S J0320-0446 (SpeX data from IBurgasser et~aTl 
I2008al) and the M8 5 dwarf HB 211 5-4518 
(jHawkins fc Besselll flQSSL ILodieu et al.l [20051 SpeX 
data from Burgasser et al. in prep.). The three spectra 
are generally equivalent in shape, exhibiting strong TiO, 
H2O and CO bands; weaker absorption from VO and 
FeH; and line absorption from K I and Na I (doublet 
lines are unresolved in these data) . The spectra of both 
SDSS J0006-0852 and 2MASS J0320-0446 also exhibit 
subtle features not present in that of HB 2115—4518, 
most notably an excess of flux at 1.27 /mi, 1.55 /im and 
(less prominently) 2.1 /im; and a distinct "no tch" feature 
at 1.6 /im. As discussed in|Burgasscr (2007Q) and shown 
in Section 4, these features indicate the presence of an 
unresolved T dwarf companion to SDSS J0006-0852. 

2.4. High Resolution Near-Infrared Spectroscopy 

High resolution near-infrared spectra of 
SDSS J0006-0852 and LP 704-48 were obtained 
on seven epochs between 2010 August and 2012 January 
using the NIRSPEC echelle spectrograph on the Keck 
II telescope. Observations are summarized in Table G3 
On four epochs we used the N3 order-sorting filter and 
0'.'432xl2" slit to obtain 1.18-1.30 /tm spectra over 
orders 56-66 with resolution A/AA = 20,000 {Av = 
15 km s -1 ) and dispersion of 0.181 A pixel -1 . On 2011 



2.5 




HB 2115-4518 (M8.5) . 






SDSS JOOOfi-0852 (M9) . 




No I 

1 Kl 


2MASSJ032O-0446(M8.5+T5:) _ 


2.0 


I " 
FeH A*jl -1 A"\ 


[[,0 




Jlj ' FeH 


B FeH H.O 




HO/ 




1.5 


VO| — 

" r J ft j\ 






Tio 1/ j yufi/ 






r y i /* 




1.0 
















0.5 






0.0 


r . - 1- .- . - .. 


. , ...1 .., .., ...1 .., ., ^ 



1.0 1.5 2.0 

Wavelength (urn) 

Fig. 4.— SpeX prism spectrum of SDSS J0006-0852 (cen- 
ter) compared to thos e of the M8.5 HB 2115-4518 (top; 
IHawkin s & Bcsscll 1988; Burgasser et al., in prep.) and 
the M8.5+T5 sp e ctral binary 2MASS J0320-0446 (bottom; 
Wilson et al. 2003; Burgasser et al. 2008a). All spectra are nor- 
malized in the 1.2-1.3 (im range and offset by constants. The 
HB 2115—4518 spectrum (in red) is also overlaid on the other two 
spectra for comparison. Major molecular and atomic absorption 
features are labeled, as is the "notch" feature at 1.6 fim indicat- 
ing the presence of an unresolved T dwarf companion. Regions of 
strong telluric absorption are indicated by the horizontal lines. 

September 7 and 2012 January 6 (UT) we observed 
SDSS J0006-0852 and LP 704-48 with the N7 filter 
and 0'.'432xl2" slit to obtain 2.00-2.39 /im spectra over 
orders 32-38 with A/AA = 20,000 {Av = 15 km s -1 ) 
and dispersion of 0.315 A pixel -1 . On 2012 January 15 
(UT) we observed SDSS J0006-0852 with the N5 filter 
and 0'.'288x24" slit to obtain 1.43-1.70 /tm spectra over 
orders 45-53 with A/AA = 30,000 {Av = 10 km s -1 ). 
Echelle and cross-dispersion gratings were set to the 
values listed in Table [3l which varied slightly from run 
to run to maintain consistent projection of arclamp 
images on the detector. We also observed seven L 
dwarfs with previ ously measured radial velocities from 
iBlake et"ail (|2010D to serve as radial velocity standards 
(Table ??). For each source, we obtained spectra in AB 
or ABBA nodding sequences. We also observed a nearby 
AO V star for flux calibration and telluric correction, 
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with the exception of some N5 observations for which 
a comparison star was not required (see below). Dark, 
quartz lamp and NeArXeKr arc lamp frames were 
obtained at the beginning or end of each night without 
changing the instrument configuration. 

All data were reduced using REDSPEC, an IDL-based 
software package developed at UCLA for NIRSPEC by 
S. Kim, L. Prato, and I. McLeanQ. Images were first 
corrected for pixel-response variations using the dark- 
subtracted flat field frames. Individual orders were then 
isolated, rectified, and pair-wise subtracted, and spectra 
were extracted by summing across rows. In this study, 
we focused exclusively on order 59 (1.283-1.300 /xm) in 
the N3 data, order 49 (1.545-1.567 fjaa.) in the N5 data, 
and order 33 (2.291-2.326 /xm) in the N7 data. The lat- 
ter two o rders are relatively devoid of telluric abs orption 
features (|Prato et all 120021: IMcLean et all 120071 ) . while 
order 33 samples the strong CO band around 2.3 /mi 
(|Blake et al.l 120081 ) . Image rectification and wavelength 
calibration were performed using telluric OH emission 
lines pre sent in the long scienc e exposures for orders 49 
and 59 (jRousselot et al.l [20001) ; for order 33, we used 
the arc lamp images and a vacuum line list from the 
National Institute of Standards and Technology ( NIST ) 
atomic spectral line database (Ralchen ko et al.l 1201 ID . 
For each observation, wavelength solutions were cor- 
rected for barycentric motion. Flux calibration and tel- 
luric absorption corrections (except for order 49) were 
calculated from the AO V spectra assuming a 9480 K 
blackbody. No filtering of fringing was performed. 



- NIRSPEC order 59 



SDSS J0006-0852 



l.i t 2010 Nov 26 

z 2010 Aug 19 

2MASS J0523-1403 



1.285 



1.290 1.295 
Wavelength (microns) 



| 1.0 



NIRSPEC order 49 

SDSS J0006-0852 



O.Or- 

1.545 



1.550 



1.555 1.560 
Wavelength (microns) 



1.565 



Fig. 6.— NIRSPEC order 49 (N5) spectrum of SDSS J0006-0852 
obtained on 2012 Jan 15 (UT, top) compared t o that of the M9 ra - 
dial velocity standard LHS 2065 (bottom; from lPrato et al.ll2002T) . 
Spectra have been corrected for barycentric motion, normalized to 
the median flux and offset by integer constants. Features arise 
primarily from H2O and FeH absorption. 



I 2 



NIRSPEC order 33 

SDSS J0006-0852 




: 2MASS J1731+2721 20111™, 10 : 



Fig. 5.— Four epochs of NIRSPEC order 59 (N3) spectra of 
SDSS J0006— 0852 (top) compared to the L2.5 radial velocity stan- 
dard 2MASS J0523-1403 (bottom). Spectra have been corrected 
for barycentric motion, normalized to the median flux and offset 
by integer constants. The 1.290 fim Mn I line is labelled; other 
features arise primarily from FeH and H2O absorption. 

Figures [5] through [7] display all of the reduced spectral 
data for SDSS J0006-0852 and LP 704-48 compared to 
select radial velocity standards. Signal-to- noise (S/N) for 
these data are generally >50, with the exception of our 
N5 data which has S/N w 20. All of the spectra show 
numerous molecular transitions arising from H 2 and 
FeH in orders 59 and 49, and CO in order 33. We also 
detect the weak Mn I line at 1.2903 am in the spectrum 
of SDSS J0006-0852 (|McLean et al.il200l . 

3. LP 704-48 AND SDSS J0006-0852 
9 See http : //www2 . keck .hawaii . edu/ inst/nirspec/redspec/index 



2.295 2.300 2.305 2.310 2.315 2.320 2.325 
Wavelength (microns) 

Fig. 7.— NIRSPEC order 33 (N7) spectra of SDSS J0006-0852 
(top two spectra), LP 704-48 (second from the bottom) and the L0 
radial velocity standard 2MASS J1731+2721 (bottom). Spectra 
have been corrected for barycentric motion, normalized to the me- 
dian flux and offset by integer constants. Features arise primarily 
from CO absorption. 



3.1. Evidence for Common Proper Motion and Distance 

Our Nickel astrometry (Figure [5]) shows convinc- 
ingly that LP 704-48 and SDSS J0006-0852 com- 
prise a common proper motion pair. To quan- 
tify this, we combined our measurements with 
prior astrometry f r om th e SuperCosmos Sky Sur- 
vey (IHamblv et al.l l2T)0ll blial) . the 2MASS All-Sky 
Point Source Catalog , and SDSS Data Release 7 
(|Abazaiian et al.l [2009). Linear fits yield proper mo- 
tions of [/j, a cos 5, [is] = [— 70±3,— 314±2] mas yr _1 
and [-84±28,-337±17] mas yr" 1 for LP 704-48 and 
SDSS J0006— 0852, respectively. Th e form er agrees with 
values reported by iSalim fc Gould! (|2003fl . These mo- 
tions are consistent with each other to within 0.5c in 
Right Ascension and 1.3c in declination. The combined 
astrometry also yield angular separations of 25'.'26±0'.'09 
in Right Ascension and — 10'.'64±0'.'17 in declination, 
pointing from LP 704-48 to SDSS J0006-0852. 

We also find that the estimated distances to these 
sources are in formal agreement, based on a com- 
bin ation of several optical and near-infrared abso- 
Mlfte magnitude/color and absolute magnitude/spectral 
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TABLE 3 

NIRSPEC Observations of SDSS J0006-0852 and LP 704-48 



Source 



Date (UT) 



MJD 



Filter A (fim) 



Slit 



Echelle/Grating (°) T (s) Airmass Calibrator 



Conditi 



SDSS J0006- 


-0852 


2010 Aug 19 


55437.5 


N3 


1.18- 


-1.30 


/ .'432xl2 


62.95/34.10 


3600 


1.13- 


-1.20 


HD 3604 


clear, seeing « 


SDSS J0006- 


-0852 


2010 Nov 26 


55526.3 


N3 


1.18- 


-1.30 


/ .'432xl2 


62.95/34.10 


4000 


1.14- 


-1.33 


HD 3604 


high humidity, f 


SDSS J0006- 


-0852 


2011 Jul 06 


55748.6 


N3 


1.18- 


-1.30 


/ .'432xl2 


62.95/34.08 


1500 


1.14- 


-1.16 


HD 3604 


clear, seeing « i 


SDSS J0006- 


-0852 


2011 Sep 07 


55811.5 


N7 


2.00- 


-2.39 


/ .'432xl2 


63.00/35.46 


2000 


1.25- 


-1.37 


HD 13936 


light cirrus, seei 


SDSS J0006- 


-0852 


2011 Sep 10 


55814.5 


N3 


1.18- 


-1.30 


/ .'432xl2 


62.95/34.08 


2000 


1.24- 


-1.41 


HD 1154 


clear, seeing « 


SDSS J0006- 


-0852 


2012 Jan 6 


55932.2 


N7 


2.00- 


-2.39 


/ .'432xl2 


63.01/35.47 


2000 


1.24- 


-1.35 


HD 1154 


clear, seeing « i 


LP 704-48 




2012 Jan 6 


55932.2 


N7 


2.00- 


-2.39 


0'.'432xl2 


63.01/35.47 


800 


1.19- 


-1.20 


HD 1154 


clear, seeing « i 


SDSS J0006- 


-0852 


2012 Jan 15 


55941.2 


N5 


1.43- 


-1.70 


0?288x24 


63.04/36.30 


2400 


1.28- 


-1.40 


48 Cet 


light clouds, see 



TABLE 4 

NIRSPEC Observations of Radial Velocity Calibrators 



Source 


Optical 


Obs. Date 


Filter 


Integration 


Airmass 


Calibrator 


RV a 


Refs 




SpT 


(UT) 




w 






(km s _1 ) 




2MASS J05233822-1403022 


L2.5 


2009 Oct 31 


N3 


2000 


1.20-1.22 


r Lep 


12.21±0.09 


1 


2MASS J09211410-2104446 


LI. 5 


2011 Oct 26 


N3 


1200 


1.32-1.33 


HD 82724 


80.53±0.11 


2 


2MASS J10220489+0200477 


LI 


2011 Mar 18 


N7 


1800 


1.28 


39 Uma 


19.29±0.11 


2 


2MASS J15150083+4847416 


L6 


2011 Sep 10 


N7 


350 


1.57-1.63 


HD 143187 


-29.97±0.14 


3 


DENIS-P J170548. 38-051645.7 


LI 


2011 Aug 11 


N7 


1200 


1.24 


HD 159415 


12.19±0.11 


4 


2MASS J17312974+2721233 


L0 


2011 Jun 10 


N7 


1000 


1.03 


HD 165029 


-29.76+0.11 


2 


2MASS J22244381-0158521 


L4.5 


2011 Sep 07 


N7 


1500 


1.08-1.11 


HD 198070 


-37.55+0.09 


5 



Ref erences. — fl UCruz et all II2003D : f2) IReid et~aT1 il2008T) : (3) IWilson et alj (I2003T) : f4) IKendall et alj d2004D : ( 51 IKirkpatrick et all 
1120001) . 

a Radial velocity measurements from [Blake et aT~ (2010). 

corporate uncertainties in photometry, spectral classifi- 

„ TABL ?iL„ „or„ cations (0.5 subtypes), and the reported systematic un- 

Distance Estimates for SDSS J0006-0852 and • . • .i L i 1 j w ct 1 

LP 704-48 certainty m the absolute magnitude relations, borne rela- 

tionships show better agreement between the two compo- 
nents than others, with the absolute magnitude/spectral 
type relations diverging the most. Using only the color 
relations, we find mean distances of 22±5 pc for LP 704- 
48 and 32±2 pc for SDSS J0006-0852, a 2cr difference. 
Note that nearly all of the relations place LP 704-48 
roughly 50% closer, suggesting that it to may be an 
unresolved binary, although age/surface gravity effects 
cannot be ruled out (see below). We adopt a weighted 
mean of 30±3 pc for the combined system, implying a 
projected separation of 820±120 AU. 

3.2. Probability of Association 

To assess the probability of chance alignment for 
thi s wide pairi n g, we followed the method described 
in iDhital et alj (|2010f) . which estimates the frequency 
of unrelated pairings using a Gala ctic model based 
on an empirical lumin osity function (jJuric et al.l 120081 : 
IBochanski et al.l I2010D and an empir ical space velocity 
distribution (jBochanski et al.1 l2007al) . The number of 
single (and hence unrelated) stars within a 6D ellipsoid 
defined by the angular separation of the binary, the es- 
timated distance to the binary, the space motions of the 
binary, and the uncertainties in these values was deter- 
mined through Monte Carlo simulation. From 10 6 sim- 
ulations, we found that only 0.0015 stars were spatially 
coincident and had proper motions similar to the values 
observed for this wide pair. When the (far more precise) 
radial motions of SDSS J0006-0852 and LP 704-48 were 



Method 




Distance (pc) 




LP 704-48 


SDSS J0006-0852 


A/a 


H02 M, vs i - z 


28±5 


31±6 


-0.4 


H02 Mi vs i- J 


31±6 


30±6 


+0.1 


H02 M T vs SpT a 


17±4 


36±7 


-2.4 


C03 Mj vs SpT a 


18±3 


35±4 


-3.4 


W05 Mi vs i - z 


17±2 


33±5 


-3.0 


W05 Mi vs i - J 


27±3 


31±3 


-0.3 


B10 M r vs r - z 


22±4 


b 




B10 M r vs r - i 


21±4 


b 




B10 Mr vsi-z 


26±6 


b 




S10 Mi vs i - z 


b 


31±8 




S10 Mivsi-J 


b 


36±6 




D12 Mj vs SpT a 


20±5 


37±8 


-1.8 


D12 M H vs SpT a 


20±4 


40±8 


-2.2 


D12 M K . vs SpT a 


20±4 


42±9 


-2.2 


Weighted Mean 


22±5 


32±2 


-2.0 



Ref erences. — ( H02): IHawlev et al]j2002D ; ( C03): ICruz et a 
1 120031) ; fW 05l: IWest et alj 1120051) ; (B1 0): IBochanski et~ 
1 120101) : (S10): ISchmidt eraT] 120101) ; (D12): IDupuv fe Lrull2012al 

a Excluded from final average. 
Outside defined color range. 



type relati o nshipjrci defined in Hawlev et al.l (I2002D; 
, et all (l2003h~IWest et~all (120051): ISchmidt et al 




IBochanski et al.1 ()2010l ): and IDupuv & Lhi 
Table [5] summarizes these estimates, which in- 

10 Note that l i near coefficient of M r versus r — i in Table 4 of 
BochanskFeOD pOlTJ) should be +4.548, not -4.548. 



7 



also considered (see below), the number of chance align- 
ments fell to zero (i.e. probability <10 -6 ). We there- 
fore conclude with high confidence that LP 704-48 and 
SDSS J0006— 0852 are not a chance alignment. 

3.3. Activity, Age and Metallicity 

Additional constraints on the physical properties of 
this system can be inferred from our optical spec- 
troscopy of LP 704-48 and SDSS J0006-0852. The 
lack of Li I absorption in either source rules out masses 
below 0.06 Mq, implying ages older than 90 Myr 
and 210 Myr, respectiv e ly, ba sed on the evolutionary 
models of iBaraffe efail |2003) and T e ff estimates of 
2660±150 K and 240 0±160 K as i n ferred from the 
Te/z/SpT relation of IStephens eFall (|2009)H More 
stringent age constraints come from the lack of Ha emis- 
sion in either source. For SDSS J0006— 0852, one could 
attribute the lack of emission to an increasingly neu- 
tral photosphere that is decoupled from the magnetic 
field, a hypothesis used to explain the decline in both 
Ha and X-ray emission in dwarfs later than M8-M9 
(e g.-IGizis et al.ll2000bt IMohantv et~aT1l2002t IWest et al.l 
120041) . However, the inactivity of LP 704-48 is remark- 
able given the high incidence of Ha emission among M7 
dwarfs in the vicinity of the Sun (> 90% of sources with 
vertical scaleheights \Z\ < 50 pc; lHawley et al.l 119961 : 
IGizis et all [2000b; We st et all 12001 12003 . 12008?^ From 
the acti vity frequencies o f M dwarfs at various scale- 
heights, IWest et all (|2008l ) have inferred an activity life- 
time for M7 dwarfs of 8t°{ 5 Gyr. Since emission from 
LP 704-48 is seen in neither our spectrum nor that of 
ICruz fc Reidl (f2002h . and given that the limit on Ha lu- 
minosity is nearly two orders of magnitude below the 
local active mean (log 10 Lual^boi = —4.3; We st et al.l 
l200l . we conclude that this source is truly inactive and 
that the LP 704-48/SDSS J0006-0852 system is likely 
to be quite old. 

This conclusion is supported by the kinematics of the 
system. Combining our proper motion measurement, 
distance estimate and systemic radial velocity determi- 
nation below, we derive the heliocentric UVW velocities 
listed in Table El These velociti es are on the b oundary of 
the young -old disk as defined in lLegge tt (1992). The disk 
classification is consistent with what appear to be near- 
solar metallcities for LP 704-48 and SDSS J0006-0852, 
as evident from their spectral energy distributions and as 
quantified by the £ index of ILepine et al.l (|2007D . Both 
sources have £ ~ 1.00 (Table (2) indicating roughly solar 
metallicities. 

We therefore conclude that LP 704-48 and 
SDSS J0006— 0852 comprise a physically associated, 
widely-separated system of relatively old, VLM stars 
with near-solar metallicity. 

4. THE M DWARF/T DWARF SPECTRAL BINARY 
SDSS J0006-0852 

As discussed in Section 2.1, the peculiar features ob- 
served in the low-resolution near-infrared spectrum of 
SDSS J0006— 0852 indicate the presence of an unseen 
T dwarf companion. To test this hypothesis, we per- 
formed a spectral template fitting analysis similar to 

11 Uncertainties include classification of 0.5 subtypes for 
both components an d a 100 K systematic uncertainty in the 
IStephens et al.l 120091) relation. 



that described in [Burgasscr et all (|2010l ). We drew 638 
spectra of 614 sources with S/N > 20 from the SpeX 
Prism Spectral Librarie^3 with published optical and/or 
near-infrared spectral types between M7 and T8, ex- 
cluding known binaries, subdwarfs and low-surface grav- 
ity brown dwarfs. The spectra were reclassified in the 
near-infrar e d usin g the index-based scheme defined in 
IBurgasseH (]2007aD . Spectral fluxes were then scaled 
to absolu te F\ unit s using the Mj /spectral type re- 
lation of ICruz et al.l (|2003f ) for types M7-L2 and the 
Mk s /spectral type relation of lLooper et al.l (|2008al) for 
types L2-T8. We then combined pairs of flux-calibrated 
spectra to create binary spectral templates (B(X)), con- 
straining the spectral types of the primary to M7-L8 
and of the secondary to L9-T8, thereby generating 44605 
unique binary templates. We compared both the origi- 
nal 638 single spectra and the binary templates to the 
spectrum of SDSS J0006-0852 (5(A)) over the wave- 
length ranges {A} = 0.95-1.35 /zm, 1.45-1.80 /xm and 
2.00-2.35 /xm using the x 2 statistic, 

^E(^P) 2 a) 

where <r(A) corresponds to the uncertainty spectrum of 
SDSS J0006-0852 alone. 

Figure [8] displays the four best-fitting binary tem- 
plates (the best-fitting single template, HB 2115—4518, 
is shown in Figure SJ . A combination of a late- type M 
dwarf and mid-type T dwarf accurately reproduces the 
spectrum of SDSS J0006— 0852, including the excess flux 
at 1.27 /an, 1.55 /jm and 2.1 /xm. This combination also 
reproduces with reasonable fidelity the notch feature at 
1.62 fim, arising from FeH absorption in the M dwarf 
primary and the pseudo-continuum peak blue ward of the 
1.6 Mm CH4 band in the T dwarf secondary (Burgasscr 
l2007bf ). Several dozen binary templates were found to 
provide statistically superior matches to the spectrum of 
SDSS J0006-0852 as compared to HB 2115-4518 alone 
(x 2 = 12.5), with significance val ues >99% based on th e 
F-test statistic (see Eqns. 2-6 in Burgasscr et al.l[2010() . 
Table [7] details the ten best binary template fits, includ- 
ing component names, near-infrared spectral types and 
relative JHK magnitudes on the MKCFl photometric 
system. Averaging over all fits with an F-test statistical 
weighting, we infer mean component types of M8.5±0.5 
and T5±l, both rounded off to the nearest half-subclass. 
The near-infrared classification of the primary is consis- 
tent with the combined-light optical classification. We 
also estimate relative magnitudes of 3.2±0.3, 3.8±0.5 
and 4.4±0.6 in the J-, if- and X-bands. The com- 
ponent properties of SDSS J0006— 0852 are essentially 
identical to thos e inferre d for 2MASS J0320-0446AB by 
iBurgasser et al.l (|2008a). which is not unexpected given 

12 See htt p : //www .browndwarf s . org/spexprism D ata we re 
draw n from Burgasscr et al. (2004, 2000a b, 2007al. l2008allbl 
[20TrT); [Cruz et all 112004); I Burgasser & McElw aml 112009); 
IChiu et ail 112009); [McElwain & Burgasscr 1 2006):~ IReid et al.l 
OOOrl); IBurgasser] f; 2007a b c); Licbcrt & Burgas sen [121)071); 
lLooper et al.l H2007I. I2008bl); JLuhman et al.l 11200 7); Sicgl er et at] 
1120071) ; ISheppard L Cushinel 1120091 ); IKirkpatrick et alj 1120101 )" 
Cruz et al. in prep.; and Burgasser et al. in prep. 

J Mauna Kea Observatory fi lter set; see ITokunaga et~ail (2002) 
and lSimons fc Tokunagal l|2002l ). 



TABLE 6 

Properties of the LP 704-48/SDSS J0006-0852AB System 



Parameter LP 704-48 SDSS J0006-0852AB System Ref. 



wpticai op i 


N 1 7 


MQ 




i ,z,o 


NIR SpT 




M8 5+0 5 4- T5+1 




1 


V 


17.313±0.005 


21.04±0.06 




4 


r — i 


2.373+0.006 


2.82±0.06 




4 


i — z 


1.197±0.006 


1.77±0.03 




4 


i-J 


1.78±0.02 


2.32±0.04 




4.5 


J 


11.97±0.02 


14.14±0.04 




5 


J-K s 


0.88±0.03 


1.01±0.05 




5 


log 10 L Ha /L bo i 


< -6.2 


< -5.7 




1 


Est. d (pc) 


22±5 


32±2 


30±3 


1 


^ a (mas yr _1 ) 


-70±3 


-84±28 


-70±3 


1 


(mas yr" 1 ) 


-314±2 


-337+17 


-315±2 


1 


RV (km s" 1 ) 


-15.3±0.9 


Variable 


-15.3±0.3 a 


1 


a (AU) 




0.286±0.009 


820±120 


1 


a(") 




0'.'0095±0"0010 b 


27'.'41±oyiO 


1 


U (km s- 1 ) 






31+3 


1 


V (km s- 1 ) 






-37±3 


1 


W (km s- 1 ) 






2±2 


1 


Est. Age (Gyr) 


>8 


^3-4 


£3-4 


1.2 


Est. Masses (M ) c 


0.092 


0.082-0.083, 0.049-0.064 


0.22-0.24 


1 



References. — (1) This pap er: (2llWest et all 120081) : f3 HCruz fc Reidl 120021) : (4) SDSS: 
lYork et all 12000D : (5) 2MASS: ISkrutskie et al.l 120061) 

a Based on constrained radial velocity orbit fit to SDSS J0006-0852AB; see Table llOl 

At maximum elongation. 
c For an age range of 3—10 Gyr, based on the evolutionary models of Baraffc et al. (2003). 



TABLE 7 

Best-Fitting Spectral Binary Templates for SDSS J0006— 0852. 



Primary 



SpT 



Secondary 



SpT 



X 



HB 2115-4518 M8.6 

DENIS-P J2353-0833 M8.6 

HB 2115-4518 M8.6 

HB 2115-4518 M8.6 

HB 2115-4518 M8.6 

DENIS-P J2353-0833 M8.6 

DENIS-P J2353-0833 M8.6 

HB 2115-4518 M8.6 

HB 2115-4518 M8.6 

HB 2115-4518 M8.6 



2MASS J21513839-4853542 T4.2 2.886 3.384 3.908 5.13 

SDSS J135852. 68+374711.9 T4.8 2.930 2.667 3.471 5.16 

SDSS J074149.15+235127.5 T5.4 3.096 3.880 4.457 5.31 

SDSS J083048.80+012831.1 T4.9 2.953 3.598 4.168 5.33 

2MASSI J2254188+312349 T3.9 3.097 3.455 3.923 5.50 

2MASSI J0937347+293142 T5.6 3.199 3.086 3.981 5.53 

SDSS J083048. 80+012831.1 T4.9 3.285 2.950 3.607 5.54 

SDSS J074201.41+205520.5 T5.2 3.107 3.838 4.492 5.54 

2MASS J05591914-1404488 T4.4 3.019 3.615 4.176 5.56 

2MASS J06020638+4043588 T4.6 2.986 3.637 4.200 5.57 



Average Primary M8.7+0.2 Average Secondary 



T4.8+1.1 3.2+0.3 3.8+0.5 4.4+0.6 99.8% a 



a Confidence level that the best-fit binary template provides a statistically better match than the best-fit single template (HB 
2115-4518, x 2 = 12-6) based on the F-test statistic. 



the similarity in their near-infrared spectra (Figure 0]). 

5. RADIAL VELOCITY VARIABILITY AND ORBIT 

5.1. Radial Velocity Measurements 

While our spectral template analysis provides com- 
pelling evidence that SDSS J0006-0852 harbors a T 
dwarf companion, we sought independent verification by 
searching for radial velocity variations in the NIRSPEC 
data. Velocities were determined by cross-correlating 
the SDSS J0006-0852 and LP 704-48 spectra with each 
of the radial velocity templates obse rved in the same 
order using the IDL xco rl package (jMohantv k. Basril 
120031: I West fc Basil |2009( ). We correlated the spectra 
over the wavelength ranges 1.283-1.299 (im (order 59), 
1.5455-1.5665 fim (order 49) and 2.297-2.325 (im (or- 
der 33), in each case sampling ten independent windows 
equally-spaced in logarithmic wavelength. The individ- 
ual correlations within a given order were combined using 



an outlier-resistant mean based on comparison of me- 
dian and median absolute deviation statistics (with a 
2.5cr threshhold), and velocities determined from multi- 
ple standards in a given epoch (Table ??) were combined 
using an uncertainty-weighted mean. We also measured 
the radial velocities for each of the standards in order 
to assess systematic effects. We found these measure- 
ments, with uncertainties ranging from 0.5-0.9 km s , 
were all consisten t, albeit wit h lowe r precision, with 
those reported in Bla ke et al.l (|2010D to within 1.6ct. 
While a forward-modeling approach would have likely 
produced better precision, our uncertainties are never- 
theless roughly equivalent to contemporar y NIRSPEC 
RV measurements of very cool dwa rf s (e.g.. IPrato et alj 
2002t iZapatero Osorio et al.l l2007t iBlake et al.l 120081 



20ialRice et al.H2010H Rodlcr et al. 2012) 



Table ?? lists the final radial velocity measurements for 
SDSS J0006— 0852 for each epoch. There is clear varia- 
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TABLE 8 

Radial Velocity Measurements of SDSS J0006— 0852 and LP 704-48 by Comparator (in km s _1 ) 



MJD Order J0921-2104 J0523-1403 J1515+4849 J1731+2721 J2224-0158 J1022+5825 J1705-0516 B] 

(80.53 ± 0.11) (12.21 ± 0.09) (-29.97 ± 0.14) (-29.76 ± 0.11) (-37.55 ± 0.09) (19.29 ± 0.11) (12.19 ± 0.11) 



SDSS J0006-0852 



55427.53279 


59 


-7.36 ± 0.38 


-7.08 ± 0.46 












55526.27187 


59 


-18.39 ± 0.50 


-17.99 ± 0.36 












55748.61963 


59 


-15.57 ± 0.93 


-15.4 ± 1.5 












55811.54273 


33 






-20.90 ± 0.65 


-20.22 ± 0.28 


-21.1 ± 1.1 


-22.17 ± 0.38 


-22.8 ± 1.1 


55814.53006 


59 


-20.17 ± 0.74 


-19.07 ± 0.85 












55932.23208 


33 






-22.60 ± 0.77 


-22.89 ± 0.45 


-23.3 ± 1.4 


-24.32 ± 0.51 


-25.35 ± 0.45 


55941.19343 


49 
















LP 704-48 


55932.23208 


33 






-13.84 ± 0.63 


-14.39 ± 0.67 


-13.6 ± 1.2 


-16.72 ± 0.33 


-15.10 ± 0.64 



Note. — Radial - 



the heliocentric 



cd Modified Julian Date (MJD; Juli; 



TABLE 9 

Radial Velocity Measurements for SDSS J0006— 0852 



Date (UT) 


MJD 


Order 


RV 








(km s -1 ) 


(km s _1 ) 


2010 Aug 19 


55427.53279 


59 


-7.25 


0.30 


2010 Nov 26 


55526.27187 


59 


-18.13 


0.30 


2011 Jul 06 


55748.61963 


59 


-15.51 


0.78 


2011 Sep 07 


55811.54273 


33 


-20.98 


0.92 


2011 Sep 10 


55814.53006 


59 


-19.70 


0.56 


2012 Jan 06 


55932.23208 


33 


-24.0 


1.1 


2012 Jan 15 


55941.19343 


19 


-21.6 


1.3 



Note. — Radial velocities (RVs) are given in the heliocentric 
reference frame for the observed Modified Julian Date (MJD; Ju- 
lian Date - 2400000). 



tion in the velocities, with values ranging from —24 to 
—7 km s^ 1 over the 16 months of observation. These val- 
ues straddle the radial velocity measured for LP 704-48, 
-15.6±0.4 km s _1 , which we consider to be an estimate 
of the systemic radial velocity (see below). The x 2 f° r 
these data, 1000 with 6 degrees of freedom, rules out a 
constant radial motion, indicating the presence of a grav- 
itationally perturbing companion. 

5.2. Radial Velocity Orbit: MCMC Analysis 

With seven measurements in hand, we have a mini- 
mum set necessary to constrain a radial velocity orbit for 
SDSS J0006-0852. To do this, we performed a Monte 
Carlo Markov Chain (MCMC) analysis following iFordl 
(2005). Wc defined the parameter vector: 

9 = {lnP,lni^i, v^cosw, Vesinw, 4>i,V CO m.}, (2) 

where P is the period of the orbit in years, K\ the semi- 
amplitude of the primary's orbital motion along the line 
of sight, e the eccentricity, u> the argument of periastron, 
4>i the mean anomaly in the first epoch (Ti) of data, and 
V com the system's center of mass radial velocity. Note 
that the epoch of per i astron passage is T = Ti — P0i/27r. 
As discussed in lFordl ()2005l ). this particular configuration 
of orbital parameters is found to improve convergence for 
MCMC analysis in cases of small e. The semi-amplitude, 



period and eccentricity are related to the mass function: 

,1/3 

/m sin i = TTTTTR = Ki 



(M!+M 2 )2/3 



2ttG 



( 3 ) 

where i is the inclination of the orbital plane with respect 
to the plane of the sky (0° corresponds to a face-on pro- 
jection), and Mi and M 2 are the masses of the M and 
T dwarf components of SDSS J0006— 0852, respectively. 
We performed three sets of 300 MCMC trials to maxi- 
mally explore the six-dimensional parameter space. For 
each trial, we started with an initial parameter set 9q se- 
lected from distributions derived from the observational 
data. For P and K\, we drew from Gaussian distribu- 
tions with means of 0.5 yr and 9 km s _1 and widths 0.2 yr 
and 3 km s _1 , derived from the scales over which the ra- 
dial velocities reversed in trend. For e, uj and <j>\, we drew 
from uniform distributions spanning ranges [0,0.9], [0,27r] 
and [0,27r], respectively. For V com , we drew from a Gaus- 
sian distribution centered on the observed radial velocity 
of LP 704-48 with a standard deviation of 0.6 km s _1 , 
under the assumption that the relative orbital motion 
of the pair is negligible (<0.5 km s _1 ) Throughout the 
simulation, we imposed boundary conditions on the pa- 
rameters of 0.05 < P < 10 yr, 1 < K x < 100 km s" 1 and 
e < 0.95. 

For each step in the MCMC chain, we computed a new 
vector 9' by changing one parameter 6j, drawing from a 
Gaussian distribution 



(4) 



where j3 is a vector of scale factors. Initial scale factors 
were chosen by trial and err or to optimize acceptance 
rates ([Gelman fc Rubinlll992| ). We also found it useful 
to allow these scale factors to increase by 20% (up to 
a factor of 20) every 500 chain steps after the associ- 
ated parameter was last updated. With each parameter 
set, we computed that model's radial velocities (V mo d) at 
the epochs of observation and the \ 2 deviation with the 
measured values (V b s ): 



x 2 = E 



V, 



mod,k 



Vobs,k 



&obs.k 



v c< 



V L p 



704-48 



0~LP 704-48 



(5) 
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Wavelength (\lm) Wavelength (fJ.m) 

Fig. 8. — SpeX spectrum of SDSS J0006— 0852 (black lines) compared to the four best-fitting binary spectral templates (purple lines), 
all normalized in the 1.1—1.3 fim region. Primary (red lines) and secondary (blue lines) component spectra for the spectral binaries are 
also shown, scaled according to their contribution to the combined-light spectrum (source names in each panel are listed). The uncertainty 
spectrum of SDSS J0006— 0852 is indicated as the grey line at bottom. Shaded regions at the top of the panels indicate the spectral ranges 
over which the spectrum of SDSS J0006— 0852 and the binary templates were compared, and \ values are listed. The inset boxes show a 
close up of the 1.50-1.75 fim region where the notch feature is detected. 
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Fig. 9. — Expected component and combined masses of 
SDSS J0006— 0852AB (solid black curves) as a functio n of age 
based on the evolutionary models of Baraffc ct al. (2003) and as- 
suming component T e// s of 2400 K and 1080 K for the M8.5 
primary and T5 secondary, respectively. The system mass func- 
tion, /m = M2/(Mi + M2) 2 / 3 (dashed line) is also shown, with 
values corresponding to the right axis. Hydrogen-, lithium- and 
deuterium-burning mass limits are indicated by the horizontal dot- 
ted lines. 

where the sum is over all observations. Note that we 
have explicitly included the measurement of the radial 
motion of LP 704-48 as an "observation" of the systemic 
motion for SDSS J0006-0852AB. The new parameter set 

0' replaced if 

C/(0, 1) < e^ 2 '-* 2 )/ 2 , (6) 

where J7(0, 1) is a random number drawn from a uniform 
distribution between and 1. This process was repeated 
cyclically for each parameter 100,000 times for each of 
the 300 trials. 

While these data are nominally sufficient to constrain 
the radial velocity orbit of the system, we also incor- 
porated constraints based on the physical properties of 
the components as inferred from their spectral types 
and e volutionary models. Using the Stephen s et al.l 
(2009) T e ///SpT relation, we adopted component effec- 
tive temperatures of 2400±160 K and 1080±170 K for 
SDSS J0006-0852A and B, respectively, and used these 
to estimate component mass es as a function of a ge with 
the evolutionary models of iBaraffe et al.l (|2003f ) . Fig- 
ure M plots these masses, the combined system mass and 
the mass function (Eqn. for ages of 0.1-10 Gyr. Note 
that the estimated mass of SDSS J0006— 0852A asymp- 
totes to 0.082 M beyond 2 Gyr, while the estimated 
mass of SDSS J0006— 0852B (and by extension the com- 
bined mass and mass ratio) increase monotonically from 
0.033 M Q to 0.064 M for ages 1-10 Gyr, never breach- 
ing the minimum mass for hydrogen burning. The mass 
ratio (q = M 2 /Mi) ranges from 0.41 at 1 Gyr to 0.78 at 
10 Gyr. The evolutionary estimates of /m constrain the 
allowed orbital solutions by the requirement that 

f°M blt sin i < maximum(/™ odei ) « 0.23 M Q /3 (7) 

(i.e., sini < 1). To explore the systematic effects of our 
external constraints on the orbit model, we performed 
separate MCMC trials with both Vcom an d evolution- 
ary model constraints (Case A), without the Vcom con- 
straints (Case B) and with neither the Vcom nor evolu- 
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Fig. 10. — Best-fit primary radial velocity orbit solution for 
SDSS J0006-0852AB based on our Case A MCMC fit. The top 
panel shows the full radial velocity curve with measurements (solid 
circles with error bars) as a function of observation date; the bot- 
tom panel shows the velocity curve as a function of phase. Open 
circles indicate the expected radial velocity at a given observational 
epoch based on the best- fit model (blue lines). Orbital parameters 
K\, P, e, and Vcom (horizontal dashed line) are listed, as is the 
X 2 value for the fit and for a constant velocity solution (x 2 )- 

tionary model constraints (Case C). 

5.3. Results 

Figure [10] displays the best-fit orbital solution for our 
Case A MCMC analysis. The observed radial velocities 
are well-fit to a low-eccentricity (e = 0.12), short-period 
(0.404 yr = 148 days) orbit, with the first and last obser- 
vations taken around maximum elongation. The inferred 
center-of-mass velocity, —15.7 km s _1 , is consistent with 
the measured velocity of LP 704-48, although this is not 
surprising for our Case A analysis given the Vcom con- 
straint. With a \ 2 = 2.21, this solution is an excellent 
fit to the data, and suggests that the uncertainties on 
the radial velocities may even be slightly overestimated. 
The best case solutions for our Case B and Case C anal- 
yses have similar parameters, notably yielding the same 
Vcom values without a constraint, and comparable \ 2 
values. Figure [TT] displays a representative relative astro- 
metric orbit for the best-fit Case A solution. The appar- 
ent separation of the two components based on the solu- 
tion never exceeds 10 mas, making SDSS J0006— 0852AB 
unresolvable with current AO or space-based instrumen- 
tation. 

Marginalized individual and joint probability distribu- 
tions for various orbit parameters and parameter pairs 
were derived by first computing a weight p for each so- 
lution i in the chain based on its x 2 value relative to the 
minimum \ 2 of the chain (i.e., the best-fit solution): 

p . = e -(xVmm[ {x 2 }])/2_ (8) 
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Fig. 11. — Predicted relative astro-metric orbit of the LP 704- 
48/SDSS J0006-0852AB system based on the best-fit radial ve- 
locity orbit, and assuming i = i m i„ and an arbitrary longitude of 
nodes. The primary is located at the origin. Both apparent (bot- 
tom and left axes) and projected (top and right axes) coordinates 
are shown, and epochs at which radial velocity measurements were 
obtained are labeled. 



TABLE 10 

Inferred Orbital Parameters for SDSS J0006-0852AB 



14 This pruning is equivalent to the common practice of exclud- 
ing some fraction of the initial steps in a given MCMC chain; see 
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Parameter 


Case A 


Case B 


Case C 


Minimum \ 2 


2.21 


2.38 


1.98 


Period (yr) 


0.404±0.004 


0.407±0.009 


0.407±0.007 


Period (dy) 


147.6±1.5 


148.7±3.3 


148.7±2.6 


Eccentricity 


0.10±0.07 


0.22±0.16 


0.10±0.06 


w (°) 


-0.3±2.1 


-67±4 


-0.1±2.2 


T (MJD) 


55421±3 


55421±5 


55421±4 


Ki (km s- 1 ) 


8.2±0.4 


8.5±0.6 


8.4±0.6 


Vcom (km s" 1 ) 


-15.7±0.2 


-16.4±0.7 


-15.9±0.5 


Su (M^ /3 ) 


0.202±0.010 


0.205±0.013 


0.202±0.010 


ai sin i (AU) 
a (AU) a 


0.111±0.005 


0.122±0.011 


0.112±0.018 


0.286±0.009 


0.288±0.010 




Minimum i (°) a 


61±5 


62±6 




Minimum Age (Gyr) a 


3.7±0.8 


4.0±1.0 





Note. — Case A incorporates constraints on both Vcom and /m ; Case B 
drops the Vcom constraints; Case C drops both Vcom and /jvr constraints; 
see Section 5.2. 

a Assuming a combined system mass of 0.144±0.013 Mq based on the evo- 
lutionary models of Burrows et al. (2001) and Baraffc et al. 2003. 

For computational ease, we excluded all parameter sets 
with a < 10 -4 , which left of order 10 5 solutions per 
CasJn. For individual parameters, we divided the full 
range of parameter values {Oj} into 30 bins, and for 
each bin summed the probabilities of all solutions with 9j 
in that bin. The resulting marginalized probability dis- 
tributions were normalized and then fit with Gaussians 
to determine parameter means and standard deviations. 
For parameter pairs, we performed the same analysis, 
dividing the ranges of both parameters into 30 bins and 
summing the probabilities of those solutions whose pa- 



FlG. 12. — Probability distributions of derived orbit parameters 
for SDSS J0006-0852AB based on our Case A MCMC fit. The 
leftmost panels show the marginalized probability distributions of 
(from top to bottom) period (in yr), eccentricity, primary radial 
velocity semi-amplitude (K\ in km s _1 ), center of mass velocity 

(Vcom m km s — 1 ), mass function (/m = M2 sini/(Mi+M2) 2/ ' 3 in 
1/3 

Mg' J ), and primary scmimajor axis (ai sine in AU). Gaussian fits 
to these distributions used to determine mean values and standard 
deviations are shown in red. The other panels display joint proba- 
bilities distributions between parameter pairs, which shading from 
light to dark blue tracing factors of 0.05, 0.1, 0.3, 0.5, 0.7 and 0.9 
relative to the best-fit solution (see Eqn.[8]|. 

rameters simultaneously satisfy both bin ranges. 

As summarized in Table 1101 we derive fairly 
stringent (<10%) constraints on the orbital pe- 
riod (0.404±0.004 yr, 147.6±1.5 days), eccentricity 
(0.10±0.07), and Ki (8.2±0.4 km s" 1 ), and determine 
oj to within 2° and periapse passage Tq to within 3 days 
(periapse occurred 6 days prior to our first observation) . 

These parameters yield f M suxi = 0.202±0.010 Mq /3 , 
which is near the maximum Jm allowed from the evo- 
lutionary models (Eqn [7]) . Our determination of Vcom 
(— 15.7±0.2 km s _1 ) is consistent with, and more accu- 
rate than, the measured radial velocity of LP 704-48. 
Inferred parameters for our Case B and Case C analy- 
ses were identical to within the uncertainties. We also 
derived the primary semi-major axis of the system 



a\ sin! 



271-^1^ 



= 0.111±0.005 AU 



(9) 



for example Ford (2005). 



for Case A, with the other Cases giving statisti- 
cally equivalent results. For comparison, the orbit of 
2MASS J0320-0446AB, which is also nearly circular, has 
a period 1.7 times longer and a primary sem i- major axis 
1.4 ti mes wider than SDSS J0006-0852AB (jBlake et al.l 
l20"lCh . Our marginalized parameter pair distributions 
show that P, e and K\ values are uncorrelated, while 
/m sini and a\ sini values show an expected correlation 
with K\. 

For our analyses that included evolutionary mass con- 
straints, Cases A and B, we determined model-dependent 
minimum inclinations and system ages fFigure[T3)). From 
Eqn [71 we derived i OT , n = 61°±5° (Case A), with so- 
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Fig. 13. — Marginalized probability distributions of the mini- 
mum orbital inclination (top, in degrees) and minimum system 
age (bottom, in Gyr) for SDSS J0006-0852AB based on our Case 
A MCMC fit. These assume component masses inferred from the 
spectral types and evolutionary models of[Baraffc et al. ( 20Q3) (see 
Figure [2}. 

lutions up to 70° being significant. However, eclipses 
remain unlikely for this pair, requiring an inclination 
i > 89? 8 [3 Constraints on system age is also based on 
the requirement that sini < 1, which implies a mini- 
mum value for the time-dependent /m (Figure [9]). Our 
marginalized probability distribution for the minimum 
age for Case A peaks at 3.7±0.8 Gyr, consistent with the 
q Gyr minimum age inferred for LP 704-48 based on 
its lack of Ha emission. Formally adopting a system age 
of 3-10 Gyr limits the mass of the secondary to 0.049- 
0.064 M , the total system mass to 0.131-0.147 M 0! 
and the mass ratio of the components to 0.60 < q < 
0.780- The total system mass can be combined with 
the period to infer a (model-dependent) semimajor axis 
of 0.286±0.009 AU. This makes SDSS J0006-0852AB 
the third tightest VLM binary discovere d to date, behind 
the s pectroscopic binary PPL 15AB feasri fe Martini 
1999; a = 0.03 AU) and the 1 Myr eclipsing brown 
dwarf binary 2MA SS J05352184-0546 085 (hereafter 
2MASS J0535-0546: IStassun et aD 12555 a = 0.04 AU), 
and tighter than 2MASS J0320-0446 (a = 0.4 AU; 
IBlake et alJ (|2"5T5h ). 

6. DISCUSSION 

6.1. Very Low Mass Triples 

The LP 704-48/SDSS J0006-0852AB system joins a 
growing list of VLM triples in which all three compo- 
nents have masses near or below 0.1 M (assuming a 

15 The geometric probability that one component will eclipse is 
Rl +R 2 l+2ej* ^ Q 4%j assuming R 1= R 2= R Jupit 

16 This is based on t he evolutionary models of IBaraffe et al.l 
12003ft . The models of ISaumon fe Marlev) <200lD yield similar 
values; the models of IBurrows et al.l 1120011 ) yield M 2 = 0.055- 
0.072 M , Mi+M 2 = 0.139-0.157 M , and q = 0.66-0.85. 



mass of 0.093 M Q for LP 704-48 based on the evolu- 
tionary models used above). Formally defining the cate- 
gory "VLM triple" as a bound system of three hydrogen- 
rich objects with a total mass of 0.3 M or less, we 
list the seven candidate and confirmed systems that 
fall into this category in Table [TT] With the excep- 
tion of LHS 1070ABC, these systems are strongly hier- 
archical, with a ratio of outer to inner separations rang- 
ing from -5 for LP 714-37ABC to ~2900 for LP 704- 
48/SDSS J0006-0852AB. They also have estimated com- 
ponent masses that are comparable to each other, al- 
though this may be a selection bias (sensitivity lim- 
its) and subject to the typically large uncertainties in 
mass estimates for substellar dwarfs. Three of these sys- 
tems have total estimated masses M to t S 0-15 M , al- 
though their tertiary components have not been solidly 
confirmed. DENIS-P J020529. 0-115925ABC (here- 
after DENIS J020 5-1159ABfCJ: iDelfosse et alJ 119971 : 
iKoerner et al.l 119991: IBouv et al J 12005^ was resolved into 
a ^7 AU binary by AO and HST im aging. Evid e nce fo r 
a third component was reported by IBouv et al.l (|2005ft . 
but only on the basis of persistent residuals in point 
spread function fitting of the resolved pair in multi- 
epoch HST images. 2MASS J08503 593+1057156ABC 
(here after 2MASS J08 50+1057ABfC); IKirkpatrick"et~aII 
fl999tlReid et~aUl2551 was resolved as a 0'.'16 binary with 
HST, and subsequently identified as a candidate triple 
based on the later classification of its primary, despite 
this component b eing roughly a magnitu de brighter in 
the near-infrared (jBurgasser etalll2011ah. Th is result 
has been called into question by IDupuv fc Liu! (|2012b[) . 
however, based on reana lysis of the combi n ed-light spec- 
trum . Kelu 1AB(C) (iRujz et all 119971: iMartin et alJ 
[1991 iLiu fe Leggettl 120051: iGelino et al.l 12006ft was re- 
solved as a 6 AU visual double with AO and HST 
imaging, and its primary identi fied as an unresolve d 
(<4 AU) L/T spectral binary by iStumpf et alJ (|2008ft , 
although there has been no evidence of RV variability 
in the combined-light system at the 1- 3 km s _1 level 
(|Basri fe Reinersll2006l: IBlake et al.ll2010ft . Given the un- 
certain nature of these triple candidates, we argue that 
LP 704-48/SDSS J0006-0852AB is currently the lowest- 
mass triple verified through multiple techniques. It is 
also the only VLM triple system in Table [TTj for which 
multiple orbits of the tight inner binary have been ob- 
served. 

The LP 704-48/SDSS J0006-0852AB system 
shares much in common with two recently un- 
covered, but slightly more massive, low-mass 
triples containing transiting subs tellar comp o nents : 
NLTT 41135AB/NL TT 41136 (llrwin et al.l [2010h 
and LHS 6343ABC ([Johnson et alJ 12011ft . Both are 
similarly composed of relatively wide M-dwarf pairs 
(55 AU and 20 AU, respectively) with one component 
hosting a tightly-orbiting (0.02 AU and 0.08 AU) 
substellar mass (0.03 M and 0.06 M ) companion. 
For these systems, the component separations are 
roughly an order of magnitude smaller than for LP 704- 
48/SDSS J0006-0852AB, but the relative inner to outer 
separations and component masses arc comparable. 
Moreover, based on evolutionary models, the tertiaries 
of these systems are also likely to be T dwarfs. The 
fact that three such systems have been identified in 
the span of three years suggests that such low-mass 
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TABLE 11 

Confirmed and Candidate Very Low Mass Triples (M to t J; 0.3 Mq) 





Components' 1 




Spectral Types 


Estimated Masses (Mq) 


Separations 


(AU) 


Ref 


A 




BC 


ABC 


ABC Total 


A-BC 


BC 





LP 714-37A 


LP 714-37BC 




M5.5 


M8 


M8.5 


0.11 


0.09 


0.08 


0.28 


36±5 


6.8±0.9 


1 


LHS 1070A 


LHS 1070BC 




M5.5 


M8.5 


M9 


0.115 


0.08 


0.077 


0.27 


~12 b 


3.57±0.07 


2.3 


LP 213-68 


LP 213-67AB 




M8 


M8 


L0 


0.09 


0.09 


0.084 


0.27 


340±60 


2.9±0.6 


1,5 


LP 704-48 


SDSS J0006 


0852AB 


M7 


M8.5 


T5 


0.092 


0.083 


0.056 


0.23 


820±120 


0.286±0.009 


6 


DENIS J0205-1159A 


DENIS J0205- 


1159B(C) C 


L5 


L8 


TO: 


0.06 


0.05 


0.04 


0.15 


7±1 


~1.3 


7 


2MASS J0850+1057B 


2MASS J0850 


-|-1057A(C) C 


L6 


L7 


L7: 


0.05 


0.05 


0.05 


0.15 


6.0±0.9 


<4 


8 


Kclu IB 


Kelu 1A(C) C 




L3p 


L0.5 


T7: 


0.05 


0.06 


0.02 


0.13 


6-4±?;f 


<4 


9 



References. — (1) Phan-Bao c tl^SWOft); (21ILeinert et alj J20pl)): (31 ISeifahrt et~aH 1 12003) ; (41 IGizis et""aT1 l!2000aD ; (5) IClose eTTal] (120031 1 : (6) 
This paper; (7) IBouv et al.l ( 1200511 ; TiT lBurgasser et al.1 H2011al) : (91 Istumpf et al.l ( 120081 ). 

a For this table, we refer to the outermost component as "A" and the inner binary as "BC" irrespective of mass or designation. 

Scaling from the inner semi-major axis using the period-mass ratio: (a ou t/ai n ) 3 — (P ou t / Pin) 2 X (M ou t/Mi n ), and values from Seifahrt ct 
al. (2008). 

c Candidate triples with unconfirmed third component. 

triple configurations may be quite common, p articularly 
among wide VLM pairs (e.g. JLaw et al.l [ 20 10), although 
a robust survey is needed to quantify the incidence of 
such systems. 

6.2. Stability of the LP 704-48/SDSS J0006-0852AB 

System 

The LP 704-48/SDSS J0006-0852AB system is at 
an extremum among the VLM triples listed in Ta- 
ble QT] in that it has both the widest outer separa- 
tion and the smallest (measured) inner separation in 
the sample. The outer pairing is remarkable for a 
system with M tot <0.25 M Q ; currently, only three 
other VLM field binaries are known to have projected 
separations >500 AllFI The gravitational binding en- 
ergy of LP 704-48 and SDSS J0006-0852AB, \E b \ < 
(2-3) xlO 41 erg, is low but not unprecedented. Re- 
cently, several wide multiples of comparab le total mass 
and binding energy h ave been found (e.g., Dhital et al.l 
l2010t[Law et al.l2010l) . including high-order systems such 
as NLTT 20346AB/2MASS J0850359+105716AB(C) 



> 




0.2 0.4 0.6 0.8 
Total System Mass (M sun ) 



1.0 



E h \ » 0.4 x 10 41 erg;lKirkpatrick et afl[T999|: [ReiTeTaT] 

J) 



VMm iBurgasser et al.l l2011at IFahertv et ahl I2OIID and 
G 124-62/ DENIS-P J144137.3-094559AB (\E b \ » 3 x 
10 41 erg: IMartm et al.l [19991 [2006t IBouv et al.l 12001 
ISeifahrt et al.l I2005D . both of which contain substcl- 
lar components. In part, it is the additional mass 
of the T dwarf tertiary that pushes the LP 704- 
48/SDSS J0006-0852AB syst em into a "normal" regime 
in mass/separation space (see IFahertv "eFaI|[20ll]) . This 
component also contributes to the long- term stability of 
the wide pair to extern al perturbation (Weinber g et al.l 
[19871: IDhitallt^[2TIT0l) . 

Because of the large difference in outer to inner sep- 
arations, LP 704-48/SDSS J0006-0852AB is also ex- 
ceptionally stable to internal dynamical disruption. In- 
ternal stability can be qu antified with the criteria of 
lEggleton fc Kiseleval ([1995), who examined the ratio of 
outer periapse a ou t{l — e Q ut) to inner apoapse ai n (l + ei n ) 
for triple star systems. For simplicity, we drop the ec- 



FlG. 14. — Stability ratio Y/Y m i n as a function of total system 
mass for a sample of low mass trip l e systems drawn from the Mul- 
tiple Star Catalog (Tbkovmin 1997); Fahcrt vet~aTI (120111 1: and this 
paper. The LP 704-48/SDSS J0006-0852AB system is indicated 
by the large r ed circle at left; th e eclipsing triples NL TT 41135AB + 
NLTT 41136 lllrwin et al.120101) and LHS 6343ABC Uohnson et al.l 
120111) are indicated by blue squares. 

ccntricity terms and approximate the ratio Y — a out /ai n , 
using as a out the projected separation of the wide pair- 
ing. For LP 704-48/SDSS J0006-0852AB, Y w 2900, 
a value several orders of magnitude greater the critical 
ratio 
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17 Koenigstuhl-l at 1800 AU 
2MASS J0126555-502239 at 5100 AU; 



and 2MASS J1258350H 
I2009D . 



4013083 at 



dCaballercl [20071. 
; (lArtigau et al.11 2007. 2009) 
6700 AU URadigan et al.l 



which is ~5 for this system, assuming qt n = Mg/M^ « 
0.7 for SDSS J0006-0852AB and q out = (M A + 
M B )/M LP 704-48 ~ 1.6(3 Figure [T41 compares the ra- 
tio Y/Y min for LP 704-48/SDSS J0006-0852AB and 
other low-mass triples listed in the most rec ent version 
(Apr il 2010) of the Multip le Star Catalog (jTokovininl 
fl997h . IFahertv et all (pOTl and Table HU By this met- 
ric, LP 704-48/SDSS J0006-0852AB stands out as the 
most internally stable low-mass triple known, slightly 

18 These estimates also imply a period ratio X = 
yl.5 ( _9o^i_^°' 5 w 10 5. 
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more stable than the N LTT 41135AB/NLTT 41136 sys- 
tem (|Irwin et all 120101 ) . With such a large semi-major 
axis ratio, we note that N-body simulations find the outer 
orbits of su ch systems are likely to h ave a significant 
eccentricity (jSterzik fe TokovininlfeOO^ . so the ratio of 
outer periapse to inner apoapse may be smaller than the 
semimajor axis ratio used here. 

6.3. On the Formation of 
LP 704-48/SDSS J0006-0852AB 

Given the now established existence of a handful of 
VLM triples like LP 704-48/SDSS J0006-0852AB, it is 
worth examining whether such systems and their charac- 
teristics match the predictions of star and brown dwarf 
formation models. We emphasize that current samples 
are far from statistically robust or complete, and detec- 
tion biases may be significant. 

Despite the low production rate of VLM multiples in 
general, current models do produce VLM higher-order 
systems. A commonly-ascribed mechanism is dynamical 
scat tering in the post-accretion phase of a young clus- 
ter (ISterzik fe purisenll2003tlDelgado-Donate et"al]|2004l: 
IGoodwin et all I2004[) . In their stud y of sm all N-body 
cluster dynamics. ISterzik fc Durisenl (|2003l ) found that 
11% of their simulated systems formed triple and higher- 
order multiples, and over 40% of the triples contained at 
least one brown dwarf component. Most of these sys- 
tems were widely separated and hierarchical, typically 
with the distant companion being of low mass (in sys- 
tems containing two brown dwarfs, the distant compan- 
ion was a brown dwarf in 90% of cases). However, all 
wide triples (a > 100 AU) with brown dwarf components 
had extremely low mass ratios (q < 0.3), whereas the sys- 
tems listed in TablefTTIare largely composed of near-equal 
mass components. Moreover, this simulation proved un- 
able to generate VLM binaries tighter than ~1 AU, and 
the authors conclude that dynamical interactions alone 
cannot explain the existence of such systems (which at 
the time included only PP1 15). Hence, it would appear 
that dynamics alone is incapable of producing VLM mul- 
tiples like LP 704-48/SDSS J0006-0852AB, a conclusion 
that has been reached for more massive higher-order mul- 
tiples as well (|Tokovininll2008D . 

Accretion plays an important role in the formation 
of multiple systems, as a mechanism for angular mo- 
mentum exchange and dissipa tion, and drivin g bina- 
ries to near-equal masses (e.g., iBate et all 120021 ) . iBatd 
(2012) recently examined VLM multiplicity in radiative 
hydrodynamic simulations of molecular cloud collapse, 
which incorporates both gas accretion (with some radia- 
tive feedback) and dynamics, albeit over shorter periods 
than pure N-body simulations due to computational con- 
straints. The Bate study produced three low-mass triple 
systems (M tot < 0.6 AU), all of which were hierarchi- 
cal with either the high mass or middle mass component 
at wide separation. While an exact clone to LP 704- 
48/SDSS J0006-0852AB was not generated, one system 
composed of a 0.15 M© plus 0.03 M©, 14 AU inner binary 
with a 0.07 M© wide companion at 194 AU is a reason- 
able analog, although the large inner eccentricity of this 
triple (e = 0.8) may ultimately lead to its dissolution. 
One fundamental shortcoming of this study is the ~1 AU 
separation range limit on cl ose in t eracti ons, again due to 
computational constraints. Bate (2009) examined VLM 



binary statistics at ten times higher resolution but over 
a shorter timeframe and employing slightly different gas 
physics. While they could successfully form 0.1-1.0 AU 
VLM binaries in these simulations, no analysis was made 
of higher-order multiples. 

The fragmentation of circumstellar disks around mas- 
sive stars is another propose d mechanism for forming 
VLM stars and brown dwarfs (|Stamatellos fc Whitworthl 
120081 120091) . Unfortunately, the predicted binary rate 
for VLM systems created in these environments is very 
low (~8% for ejected systems) and no low-mass triples 
have been created in simulations to date. Furthermore, 
while the separation distribution of simulated substellar 
pairs peaks in a range that is close to that inferred for 
SDSS J0006-0852AB (0.3 < a < 0.6 AU), these binaries 
tend to have very large eccentricities (e > 0.7), a fea- 
ture not seen in any of the VLM spectroscopic binaries 
identified to date. As such, disk fragmentation does not 
appear to be a viable mechanism for making triples like 
LP 704-48/SDSS J0006-0852AB. 

It is important to emphasize that current mod- 
els fall short in reproducing very short-period, low- 
eccentricity VLM pairs like SDSS J0006-0852AB and 
2MASS J0320-0446AB. Such systems are necessar- 
ily the product of dynamical and dissipative evolu- 
tion, as opacity-limited fragmentation constrains the ini- 
tial separations of self-gravitating masse s t o >10 AU 
([Larson! [1961 lLow fc Lvnden-Bell 119761 ). IBate et all 
(2002) resolves this problem through a combination of 
dynamical interactions and accretion from a circumbi- 
nary/circumtertiary disk, although again these simu- 
lations fail to produce stable VLM multiples. In- 
deed, three-body enc ounters more often tha n not dis- 
rupt VLM multiples (jMoeckel fc Bate! 12010). A more 
intriguing mechanism for SDSS J0006— 0852AB is Kozai- 
Lidov eccentricity perturbations induced by LP 704- 
48, follo wed by circ u larization throug h tidal friction 
(KC TF: IKozail U961 iHarringtonl Il96l IKiseleva et all 
119981: iFabrvckv fc Tremainell2007t [Socrates et al.ll2012j) . 
This mechanism has been proposed to explain the high 
fraction of spectroscopic binaries wit h tertiary com- 
panio ns (96% for SBs with P <3 d y; iTokovinin et~all 
[2001 see also iRaghavan et all [2010h . KCTF predicts 
end states with roughly circular inner orbits and large 
perio d ratios (X > 10 3 — 10 4 ; IFabrvckv fc Tremaind 
[20071 ). similar to SDSS J0006-0852AB. Unfortunately, 
the timescale for eccentricity pumping is long for LP 704- 
48/SDSS J0006-0852AB (r~P^ t /P jn 3 Gyr), while 
tidal circularization is extremely inefficient at the cur- 
rent orbital separation of SPSS JOOO6-O852AE0 dHutl 
H98lUMeibom fc Mathieull2005D . Thus, it appears that 
KCTF is probably not responsible for the current config- 

19 Using Eqns 50 and 52 from IHuH 1 119811 ). the timescale for 
circularization: 

O0024 f^Y«L Po 
kq \RiJ t 

where k is the primary's Love number, q the mass ratio, ao and Pq 
the final semimajor axis and period, R± the primary's radius and r 
the tidal lag with respect to the displacement between primary and 
secondary. Assuming k ~ 0.5 (i.e., comparable to Jupiter and Sat- 
urn), q = 0.7, ao/Ri = 0.29 AV/Rj up iter « 620 and P = 0.4 yr 
and r PS O.lPo, we derive t c i rc fa 6xl0 20 yr, far too long to have 
played a role in the dynamic evolution of SDSS J0006— 0852AB. 
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uration of this system, nor presumably for the isolated 
tight VLM binaries PP1 15AB, 2MASS J0535-0546AB 
or 2MASS J0320-0446AB. 

A final possibility is that LP 704-48 has played no 
role in the evolution of the SDSS J0006-0852AB bi- 
nary, which could have shrunk and circularized through 
dissipative interactions with an early circumbinary disk 
or a closer encounter resulting in an ejection. The 
wider pairing with LP 704-48 may simply be a normal 
(albeit rare) outcome of cloud fragmentation or paral- 
lel t rajectories during the dissolution of the natal clus - 
ter (jKouwenhoven et all 120101 : iMoeckel fc Bate! [2010). 
LP 704-48 and SDSS J0006-0852AB may be more akin 
to stellar cousins that stellar siblings. 



7. SUMMARY 

We have identified a hierarchical triple system, LP 704- 
48/SDSS J0006-0852AB, composed of a tight {a m = 
0.28 AU) M9 plus T5±l binary straddling the hydrogen 
burning mass limit, with a wide (a ou t ~ 820 AU), inac- 
tive M7 co-moving companion. The inner pair is found 
to be both a spectral (blended light) binary and a ra- 
dial velocity variable. By combining model-dependent 
mass constraints from the component spectral types, 
the radial velocity of LP 704-48, and orbit model fits 
of the RV data (spanning more than three orbital pe- 
riods), we have determined the orbital parameters of 
SDSS J0006-0852AB, including constraints on the or- 
bit inclination and the system's age (>3 Gyr). The 
latter is consistent with the age inferred for LP 704- 
48 based on its lack of Ha emission (8jZj*Q Gyr). 
SDSS J0006-0852AB is currently the third tightest VLM 
binary known, and its identification verifies that the spec- 
tral binary method can probe a relatively unexplored 
separation range among VLM multiples. With LP 704- 
48, this system forms the lowest- mass triple identified 
to date whose tertiary is unambiguously detected, and 
is one of seven confirmed and candidate VLM triples 
whose total mass is <0.3 M Q . While current star forma- 
tion models are capable of producing triples with sub- 
stellar components, we find that they do not yet pro- 
duce systems that replicate VLM triples like LP 704- 
48/SDSS J0006-0852AB. In particular, dynamical in- 
teractions alone do not appear to be responsible for 
the close separation of the SDSS J0006-0852AB pair, 
while disk fragmentation and hydrodynamic models have 
yet to produce comparable systems. We also rule out 
Kozai-Lidov perturbations and tidal circularization for 
the current configuration of this system. As most of the 
P < 1 yr VLM binaries found to date are isolated pairs, 
we speculate that LP 704-48 and SDSS J0006-0852AB 
formed as two separate systems bound only by proximity 
and/or common motion after cluster dissolution. Assess- 
ing whether this is the primary mode for VLM triple for- 
mation, and whether VLM triples are common or rare, 
will require a more statistically robust survey. Neverthe- 
less, these systems are benchmarks for studying multiple 
star formation theory; brown dwarf evolutionary models; 
and age, activity and metallicity diagnostics for VLM 



stars. 
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